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ABSTRACT 


A simulation study of a simplified model of a new and novel high-coulomb 
transfer switch has been performed. In contrast to the conventional triga- 
tron switch in which the currents are constricted by the z-pinch mechanism, 
the new switch operates in an inverse-pinch geometry formed by an all-metal 
chamber, greatly reducing hot-spot formations on the electrode surfaces. Ad- 
vantages of the new switch over the conventional switches are longer useful 
life, higher current capability and lower inductance, which improves the 
characteristics required for a high repetition rate switch. The simulation 
study performed here determines the design parameters by analytical computa- 
tions and comparison with the experimentally measured risetime, current han- 
dling capability, electrode damage, and hold-off voltages. Results from the 
simulation study show that the parameters of initial switch design can be 
determined for the anticipated switch performance. Results are In agreement 
with the experiment results. Although the model was simplified, it was 
accurate enough to determine the switch characteristics such as risetime, 
current handling capability, electrode damages, and hold-off voltages. 
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I. INTRODUCTION 


The development of an improved high-power, high-repetition-rate switch is 

essential for a number of applications, including excitation of pulsed high 

energy lasers , generation of energetic particle beams, magneto-plasma-dynamics 

(MPD) thrusters, inertial confinement fusion (ICF) , and electric launchers (also 

called rail guns). Although various high power closing and opening switch con- 

(1 2 ) 

cepts ’ have been studied, their performance demonstrated that the concepts 
still need improvements to overcome several significant state of the art limita- 
tions to high-power switching. These limitations include slow rise time (low 
dl/dt), current handling capability, recovery time, thermal dissipation, irre- 
versible switch damage, high-voltage hold-off, and various combinations of these. 
These are typical switch requirements found in special applications as mentioned 
above . 

The most widely used high-coulomb transfer switches are trigatron spark 
(3) 

gaps and thyratrons . Thyratrons have lower limits on voltage holdoff and 

current capability and therefore spark gaps are preferred over the limits 

of thyratron. The spark gaps suffer severe irreversible damage by hot-spot 

formations on the electrodes due to the z-pinch mechanism which constricts 

(4 5) 

the current to an area of a few millimeters in diameter * . Minor improve- 

ments of the geometry and/or adoption of special materials are insufficient to 
alleviate the effects of the z-pinch mechanism. The z-pinch mechanism results 
from the ponderomotive force F = J x B where J is the current density, and B 
the magnetic induction due to the current. (See Fig. 1). 

In a recent development in a spark-gap* switch, a laser beam was focused to 
trigger a spark between electrodes instead of using a trigger pin Even 

though a laser triggered switch may have a much faster rise time than the conven- 
tional spark-gap trigatron, the electrodes 1 surface damage by hot-spot formation 



due to the z pinch mechanism is unavoidable. The consideration of fundamentally 

different approaches or ideas is accordingly important to overcome the state of 
the* art limitations. 
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Figure 1. Z— pinch Mechanism and its switch configuration. 
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II. DESCRIPTION OF THE NEW SWITCH 


II-l. PRINCIPLE 

The new switch described here has a novel geometry in which an inverse- 

—V — K -v- 

pinch mechanism is operative and the ponderomotive force F = J x B is utilized 
to "disperse" and "move" the current sheet over the large area of the electrode 
surface instead of constricting the current. Figure 2 illustrates the inverse- 
pinch mechanism. Note that here the currents form the shape of a mushroom. This 
switch will be a "closing" switch, not primarily designed for "opening". 

Figure 3 shows the cross section of the new switch. The inner, outer and 
return electrodes are cylindrical in general and are placed coaxially. The 
switching of the current takes place in the vacuum chamber formed ,by the inner 
and outer electrodes. The dotted lines indicate the currents at the time of the 
initial breakdown, I q , at a later time (~1 micro-second), 1 ^, at the time of the 
peak current, 1^, and at the time when the current decreases near zero, 1^. The 
motion of the current sheet depicted here is the result pf the inverse-pinch 
mechanism working on the plasma through which the current flows. 

I 1-2. ADVANTAGES 

The following advantages over the conventional trigatron switches are recog- 
nizable from the description of the new switch above. (a) Reduced wear of the 
electrode material ; This is expected from the fact that the strong pinch of the 
current sheet is eliminated from the new switch since the inverse— pinch mechanism 
is utilized for moving the "feet" of the current over a large area. (b) Reduced 
wear of the insulator material : For the same reason as above, fast detachment, of 
the current sheet from the insulator shortens the period of exposure to the hot 
plasma and also alleviates a direct collision of ionized particles for the full 
discharge period. (c) Longer useful life : The advantages of (a) and (b) translate 
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Figure 3. Cross section of inverse pinch switch and schematic. 
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into a longer useful life for the new switch in comparison with the conventional 
switches* (d) Higher current capability : The advantages of (a) and (b) also re- 

sult in a higher current capability for the new switch in comparison with the 
conventional spark gaps* High power lasers applicable for welding, machining or 
laser weapons require switching of over 100 kA peak current at a repetition rate 
higher than 10 Hz. Due to the reduced current density and rapid motion of the 
"feet" of the current on the electrode surface, a significant increase of the 
current capability for the new switch is expected. (e) Reduced switch inductance : 
The inductance of a unit length of a coaxial current sheet is expressed as 

(y/2) £n (r /r . ) where y is permeability, r radius to the outer current, and r. 

o i o i 

radius of the inner current. For the conventional trigatron, the ratio r^/r^ 
becomes very large due to the reduction of r^ by z-pinch mechanism. For example, 

- 10 cm and r^ = .1 cm or r Q / r ^ = 100 is typical for the conventional triga- 
tron while r /r. < 2 is estimated for the new switch. A reduction of the induc- 
o 1 

tance by a factor of ^ 7 can be seen in this example. Since the risetime has the 
relation t^ = — / LC, a reduction of the component inductance is highly desirable 
for the fast pulsed power systems designed for laser excitation or generation of 
energetic particle beams. 

II-3. STATEMENT OF THE PROBLEM 

The purpose of this study is to contribute to the development of an inverse- 
pinch switch which is predicted to operate beyond the present state of the art 
limitations. The results of preliminary experiments were reported by Choi and 
Lee . ..The preliminary experiments, which were a proof of principle, were per- 
formed at half atmospheric pressure of N 2 gas in the switch chamber. The results 
from the tests were very promising. For instance, for a given 27 kilovolt trigger 
voltage and 10 kilovolt hold-off voltage, the current sheet formation was obvious 
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even though a limited number of trigger pins (4 pins in these tests) were used. 
The risetime measured by using a Rogowskii coil was less than or equal to 1 y 
sec as predicted. 

An effort to develop a switch beyond the state-of-the art limitations 
requires an analysis tool for better understanding of the switch performance 
and characteristics. The simulation model serves as an analysis tool which 
can be used to determine the problems and their solutions, to predict the switch 
performance, and to validate the parameters for a switch design. The perform- 
ance analysis using the simulation model is desirable for a better understand- 
ing of the switch characteristics and its comparison with experimental work. 

The switch characteristics are related to the operational parameters such as 
the 1 current distribution, the gap distance, the gas pressure, and the trigger 
mechanism for closing. 

The simulation model can also predict the switch performance character- 
istics, such as: the wear pattern of electrodes and insulators (due to hot- 

spot formation on electrodes and long-standing exposure of insulators to the hot 
plasma); the risetime; the high current capability; the high voltage hold- 
off, and the recovery time. Therefore, the simulation model would also help 
estimate the useful life time of the switch. 

The parameters, such as gap distance between electrodes, gas pressure, trig- 
ger mechanism, and current-sheet behavior, useful for the prototype switch 
design may be obtained by correlating the input parameters and output results. 
This approach is an effective and economic way to assess the feasibility of a 
new switch for high-power scaling. The design parameters can also be determined 
by experimental results. However, in most cases an experiment turns out to be 
a time-consuming and costly approach. ■ 

The current simulation model is limited to a simplified version and uses 
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the continuity, momentum, and magnetic field equations for the plasma 
motion with appropriate initial conditions (ICs) and boundary conditions (BCs) 
for the geometry of the inverse-pinch switch shown in Fig. 3. If coupled to 
the equation of the switching circuit (see Fig. 4), the model can be used for 
the analysis of major phenomena while the switch is closing. 

X- 

Among the major phenomena which appear while switching are the current 
sheet formation at a particular time and its motion. These are the key factors 
that determine directly or indirectly the wear patterns of electrodes and 
insulators, the current capability, the risetime and eventually the useful 
life of a switch. 

The simplified model serves to determine a few important parameters such 
as current sheet formation and motion during switch operation. However, it 
lacks the capability to characterize the switch operation with respect to the 
input and plasma parameters (see Fig. 5). Accordingly, for detailed under- 
standings of the switch performance, a full scale 2-dimensional simulation 
model is being developed using the alternating direction implicit (ADI) finite 
difference scheme which is a powerful numerical method. 

i The full-scale model consisting of the mathematical model and simulation 
software will serve as a tool for diagnostics and performance analysis, and 
as a method for a parameter estimation for the design of a prototype switch. 
This approach is necessary for developing an improved high-power, high rep-rate 
switch. Therefore, the efforts in developing the model should be carried out 
in parallel with the progress of the switch design. 

Further applications of the simulation model are possible since the 
coding is general and complete. For instance, with a. slight modification to 
t;he model, it can be used for the analysis of various plasma devices, such 
^s plasma light sources for laser pumping, and dense plasma production systems 
(i.e., hypocycloidal pinch (HCP) and z-pinch apparatus). 
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III. DESCRIPTION OF THE APPROACH 


The frame of the simulation software system was initiated with a complete 
set of plasma fluid model (MHD) equations. This set of equations satisfies 
given spatial and temporal domains which are continuous and bounded. When 
the domain of switching action is defined (for example. Figs. 3, - 7, and 8) the 
initial and boundary conditions can be identified by switch operation and con- 
figuration, respectively. In other words, the initial conditions are deter- 
mined by the input power, external electric circuit, the boundary conditions 
by the switch configuration, and prescribed conditions (such as flow field, 
thermal field, and magnetic field). Such equations with ICs and BCs can then 
be used as a mathematical model when certain assumptions are made because of 
geometrical symmetry, negligible quantity of a certain model property, computa- 
tional difficulty, and negligible effects of a certain physical phenomenon. 

The MHD equations were set up for the mathematical model with an axial 
symmetric condition (2-dimensional). Subsequently the mathematical model was 
written into numerical form using the alternating direction implicit (ADI) 
finite difference scheme and coded. 

A simplified version of the code was derived by considering the external 
electrical load circuit and the momentum equations. This simplified model has 
a capability to compute the motion of the current sheet due to the pondero- 
motive 'force, the velocity profiles, the current capability, and the inductance. 

An understanding of the mechanism of current-sheet formation and the 
current variation along the axial and radial directions with respect to time 
is absolutely necessary for the analysis of switch performance. The results 
of the study, together with the solution to an external circuit equation 
using time dependent plasma inductance and electrical parameters of the circuit, 
determine not only the risetime, but also the total current flowing at 
a particular time. 
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The current-sheet activities will also help determine the loci of the 
damage on electrode surface along the axial direction and thereafter the ther- 
mal dissipation from the hot-spot. 

There are three regimes of switch operation: 1) closure, or initiation 

of the switch conductive phase; 2) conduction or rundown phase; and 3) re- 
covery, or opening phase. The second (conduction) phase that occurs along the 
annular portion (the Mather geometry) becomes an important part because its 
results would show a picture of the current sheet formation and motion and the 
total current flowing as a time dependent function. This analysis is somewhat 
similar to a plasma focus problem 

The closure phase plays an important role in determining the risetime, and 
the jitter effects. The trigger pins or laser trigger form a non-uniform azi- 
muthal distribution of the initial arc or spark between triggers until plasma 
growth is fully achieved. In this case, the geometry involved in the analysis 
would be spherical coordinates (3 dimensional) as opposed to the cylindrical 
coordinates (2 dimensional) in the second (conduction) phase. Due to the geo- 
metry, non— uniformity of a plasma, jitter, and etc., the first phase for closure 
is a complicated problem. Hence, none of the analysis has been done in this 
area. 

The third phase, or opening, takes place when the current sheet is focused 
and constricted due to z-pinch mechanism on the top of the inner electrode 
shown in Fig. 3. To avoid the dense plasma focus on the top of the inner elec- 
trode, the opening phase for the purpose of a closing switch is designed to 
carry a diminishing current after the peak current. Hence, the spot damage 
on the electrode is avoided or alleviated. An opening switch version which 
uses the principle of m=o instability to cut off the current during the plasma 
focus is shown in Fig. 6. This opening switch has only a slightly different 
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configuration to the inner electrode than that shown in Fig. 3. Therefore,* this 
inverse-pinch closing switch can be easily converted into the opening switch by 
a minor design change of the electrode. One possible geometry is shown in 
Fig. 6. ' 

The analysis of the third phase will determine the recovery time. However, 
as in the first phase, the third phase poses as a complicated problem which has 
not similarity with the second phase. Hence, the" problem must be treated with 
boundary conditions which are related to the Z-pinch mechanism and coupled with 
the second phase. The results of the second phase. could be used as the initial 

■f , 

conditions to the third phase. ’ 

The following MHD equations are the expressions that were used for the 
mathematical model of the inverse-pinch high-power switch: 

The continuity equation : 

|e + v.(pv)-o . - ’ (i) 

The momentum equation : 

p + P(V.V) V +. V [n (k T "+ k T ± )]: = ^(JxB) - V-7T (2) 

The energy equation for ions : 

k 3 k ^ Vot 

7T 3F (n T i ) + w V ' (n T i V) + k nT t (V-V) + ,TT aj! ^ + V-q . 

■ Q t + § W'VTj +.pV- (i)] (3) 

The energy equation for electrons : 

Y^l lt (n T e ) + y^T V * (nT e V) + knT e (V-V)+V-cl e = ~\ + ^ j2+ ^ ei (*> 
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The magnetic field equation : 


I? = Vx(VxH) - Vx(VJ) - -f Vx ( JxH ^ VP ) 
o t 2 p 


( 5 ) 


The above equations are theoretically well defined and can be found in 
almost every plasma physics book. Therefore, further explanation about the 
terms in the equations is omitted here for brevity. 

Equations 1-5 in their applications to the switching domain are inte- 
grated with the external electric circuit equations. The electric circuit 
equations were developed based on a standard circuit for a switched power supply 
(see Fig. 4). The terminal voltage drop (V ) at capacitor is 


V 

c 


j T dl dl 

R I + L 4- L 3-2- 4- L — + R I . 

E E dt s dt p dt p p 


( 6 ) 


and the terminal voltages for the leak and plasma current during switching 
are the same. Thus 


dI i 

r l x l + l l dT 


(L + L ) 
s p 


dl 

-r- 11 + R 1 
dt p p 


(7) 


The total circuit current is 


dV 

1 = 1+1 = ~ 

L p dt 


( 8 ) 


There is another equation 


(9) 


that expresses the rate of change of the 


plasma inductance. The plasma inductance (L ) varies with respect to time, 

P 


dL Q R rz zV 

T- 2 = 2x10 9 [V In — + V -y 2 + — -] 
dt z r z „2 r 


(9) 
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where the subscripts o, r, t, z denote the outside diameter, r-direction, 
bottom tip of a node, and z-direction, respectively. 

When the above 9 equations are transformed into numerical equations by 
employing ADI scheme, the number of equations increases to 14. The solutions 
to the 14 numerical equations require the inverse of tri-diagonal matrices which 
can be done by the TRI-DIAGONAL routine more efficiently and faster than the 
conventional matrix inversion. As a whole, the ADI scheme and TRI-DIAGONAL 
routine will be the main routines for handling the model equations with reason- 
able accuracy and stability. 

The general expressions of the plasma parameters were already included in 
Eqs . 1-5. These parameters are the ion viscosity tensor, (tt in Eqs. 2 and 3), 
the conduction loss of ions (q^ in Eq. 3), the conduction loss of electrons 
(q in Eq. 4), the Hall effects (the last terms of Eqs. 3 and 5), the ion 
energy change due to electron-ion collisions (Q^ in Eqs. 3 and 4), the Joule 
heating (r| in Eq. 4), the Bremsstrahlung (W ^ in Eq. 4), the electron density 
(n in Eqs. 2, 3 and 4), the magnetic visocity (v in Eq . 5), and the Coulomb 
logarithm. 

The precise expressions of plasma behavior and switch performance are 
essential to analyze the switch characteristics, and to define and identify 
the input and output parameters for the prediction and estimation for the 
switch design. To obtain the precise expressions of plasma behavior and switch 
performance, the integration of important plasma parameters into the mathemati- 
cal model is required. However, in the simplified model, these parameters are 
ignored because of the difficulty of computation or in some cases they have a 
negligible effect. 
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Since the momentum equation can approximately describe the current formation 
and motion, the energy equations for electron and ion were omitted in the 
simplified model. Accordingly, the simplified model consists of two momentum 
equations for axial and radial flow, and four (Eqs. 6, 7, 8 and 9) for external 
electrical load circuit. 

Rewritting Eq. 2 for the radial and axial direction variations of the 
current sheet due to the ponderomotive force, we have 

pfr + P V I7 + pu|r = - J 2 B o (10) 


for the momentum equation in the radial direction, and 


+ p V I? + p V Tz • J r B 6 


(ID 


for the momentum equation in the axial direction. 

First, for defining the plasma inductance Eq. 9 is solved with the initial 
velocities. Then Eqs. 6, 7 and 8 with the plasma inductance (which was obtained 
from Eq. 9) are solved to determine the switch current, leak current, and termi- 
nal voltage drop at capacitor. The results from Eqs. 6, 7 and 8 are then used 
in Eqs. 10 and 11. In Eqs. 10 and 11 the current densities (J z and J^) are 
represented by the induced magnetic field (J = VxH) . Accordingly, the switch 
current is used to compute the induced field strength in the switch chamber. At 
each time step, loci of the current sheet front are determined by the velocities 
obtained from Eqs. 10 and 11, and at the same time the plasma inductance is deter- 
mined and used as input data in the external circuit Eqs. 6, 7 and 8. The new 
< 

switch current is then located on the front of the current sheet and used to 
determine the velocities and subsequent current sheet front. 
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IV. RESULTS AND DISCUSSIONS 


In the simplified model the imput parameters such as inductances, resis- 
tances, and capacitance of the external circuit (Fig. 4) and switch configuration 
(Fig. 3) were selected based on the experiment performed earlier The para- 

meters for external electric circuit are the followings: 


Capacitor charge 

20.xl0 ^ 

Farad 

Hold-off Voltage 

15 kV 


Switch Inductance 

41.xl0" 9 

*v 

Henry 

System Inductance 

700xl0~ 9 

Henry 

System Resistance 

1 . 0 ohm 


Leak Circuit Resistance 

5 . 0 ohm 



The switch geometry are the followings: 

Radius of Inner Electrode (Anode): 0.0365 m 

Radius of Outer Electrode (Cathode) : 0.0643 m 

Length of Inner Electrode: 0.2081 m 

Length of Outer Electrode: 0.235 m 

And others are: 

Gas Pressure: 7 mTorr 

Gas Density: 5.84x10 ^ kg/m^ 

With the above parameters the simulation procedure was set up to test the 
switch characteristics by varying the system inductance and the gas pressure. 

The system inductances selected here for the test are 100 nH, 400 nH, and 
700 nH, and the gas pressures are 0.7 mTorr, 4.3 mTorr, and 8 mTorr. The 
system inductance determines the risetime. The risetime and the current, sheet 
velocity eventually determine the loci of peak current on the electrode. Fur- 
thermore, the inductance and current flow pattern determine the exposure period 
and the wear of the electrodes and insulator. The gas pressure 'is a factot that 
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also determines a hold-off voltage, plasma motiop and traveling distance of 

plasma which will eventually determine the switch size. At relatively high 

/* 

gas pressure, the distance the current sheet trayels becomes short, so that 
the discharge takes place over a limited area of the electrodes. In the case 

of high pressure and high current flow the damage on the electrodes would be 

/■ 

inevitable even though the inverse-pinch is operative. However, the situation 
is still better when compared to the z-pinch mechanism. 

The results obtained from the simplified model by manipulating the previous 

$ 

two parameters, the system inductance and the gas pressure, are shown in Figs. 
9-32. Several parameters were also selected and are plotted which describes 
the switch performance. 

Figures 9, 10, and 11 show that the loci of the current sheet vary with 
the gas pressure in the switch chamber and with the system inductance. In 
Fig. 9 the location of the current sheet was calculated along the radial 
direction between the electrodes while the gas pressure (p) was maintained 
at 0.7 mTorr and the system inductance (L g ) at 100 nH. In Fig. 10, the 
pressure was at 8 mTorr with' the same system inductance as in Fig. 9. The 
density or pressure in the switch chamber dramatically controls the current 
sheet motion. At a system inductance of L g = 100 nH the risetime was approxi- 
mately 0.8 ]i sec as shown in Fig. 12. During the risetime the progress of the 
current sheet was limited to about 6 cm from the bottom of the switch for a 
gas pressure p.= 0.7 mTorr (see Fig. 9) and to about 2 cm from the bottom of 
the switch when the gas pressure was 8 mTorr. When the energy equation was 
included to take the shock. front into account, the real front of the current 
sheet would be slight above the location of the current sheet calculated from 
the simplified model. Figure 11 shows^-the location of the current sheet when 

the system inductance was L = 700 nH. In this case the risetime is about 

. sys . 


17 


3 p sec as shown in Fig. 13. Accordingly, the location of the peak current 
flow can be seen at about 22 cm from the bottom of the switch. However, the 
length of the actual anode is 20.81 cm so that the peak current occurs in the 
region of the z-pinch phenomena on the top of the anode electrode. 

Figures 12 and 13 show that the discharge switch current varies with the 
system inductance rather than the effects from pressure changes. The chamber 
gas pressure has a great influence on the discharge plasma inductance. How- 
ever, the effects on the discharge current from the gas pressure are not visi- 
ble in Fig. 12 because the discharge inductance is negligibly small compared' 
to the system inductance. 

Figures 14, 15 and 16 show the axial velocity variations at 5 selected 
radial positions during the discharge period. The velocity variation in 
Fig. 14 was obtained at the gas pressure p = 0.7 mTorr and the system induct- 
ance was L s y s = 100 nH. In Fig. 15, the gas pressures were 8 mTorr. The 
differences shown in Fig. 14 and 15 were due to the gas pressure. In Fig. 16 
the pressure was 0.7 mTorr and the system inductance was 700 nH. As illustra- 
ted in Figs. 14, 15, 16, 17, 18 and 19 the pattern of velocity variations re- 
sembles the discharge current variation which is determined predominantly by 
the system inductance. Figures 17, 18 and 19 show the variations of radial 
velocity during the discharge period. Figure 13 shows that by reducing the 
system inductance the risetime can be shortened by a noticeable amount. 

Figures 20, 21, 22 show the axial velocity distributions along the radial 
direction at 3 selected time. The axial and radial velocities are determined 
predominantly by the system inductance, the discharge current, and the gas 
pressure in the switch chamber. Figure 23 shows the axial velocity distri- 
butions as a function of gas pressure. Figure 24 shows the axial velocity as a 
function of the system inductance. In the same way Figs. 25 and 26 show 
the radial velocity distributions along the radial direction. 
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Figures 27 and 28 illustrate the variations of leakage current. Figure 
26 shows that the leakage current is little affected by the change of gas 
pressure. However, the system inductance has an influence on the leakage cur- 
rent as shown in Fig. 28. In this case a leakage resistance of 5 ohms was 
used. Figure 29 illustrates the discharge voltage variations with respect to 
the system inductances during the discharge period. Figure 30 shows the vari- 
ation of plasma inductance at variable gas pressure during the discharge per- 
iod. Alghough the plasma inductances are increased due to the motion of the 
current sheet, the plasma inductance is several orders of magnitude less than 
the system inductance. Hence, to achieve a fast risetime while switching, 
an effort to decrease the system inductance must first be made. 



V. CONCLUSION 


Even when a simplified model is used to simulate the inverse-pinch switch 
shown in Fig. 3, the results clearly show that the simplified simulation model 
performs well enough to describe switch performance and characteristics. The 
results shown in the figures compare well with results from the experiment . 
For instance, in the experiment the risetime with a system inductance of 700 nH 
and a gas pressure of 0.7 mTorr was observed to be about 3 y sec , which is very 
close to the result of the calculation (~ 2.75 y sec) shown in Fig. 13. 

As discussed in the section IV, the inverse-pinch switch mechanism is 
consistent with the simulation study and the results of experiment. The re- 
sults also show that the inverse-pinch switch can overcome the state of the 
art limitations, by its advantages such as reduced wear of electrodes and in- 
sulator, fast risetime, high current capability, and consequently a longer use- 
ful life. 

A new inverse-pinch switch was designed and built during the period of 
the contract. This new switch is currently installed and being tested in the 
laboratory. The detailed design of the switch is illustrated in the assembly 
and parts drawings which are attached in Appendix A, 
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C - EXTERNAL CAPACITANCE 

R_ - EXTERNAL RESISTANCE 
E 

L - EXTERNAL INDUCTANCE ■ 

E 

L - 'LEAK INDUCTANCE 
L 

R - LEAK CURRENT RESISTANCE 

Lt 

L - SWITCH INDUCTANCE 

o 

L p - PLASMA. INDUCTANCE 
I - TOTAL CIRCUIT CURRENT 
I - LEAK CURRENT 

L 

I p - PLASMA CURRENT 


Figure 4. External Switching Circuit Diagram. 
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LOAD 


INPUT PARAMETERS: 

: k 

CIRCUIT RESISTANCE (a & b) 
CIRCUIT INDUCTANCE (a & b) 

LEAK CURRENT RESISTANCE (a & b) 
LEAK CURRENT INDUCTANCE (a &. b) 
SWITCH INDUCTANCE (a, b, & c) 
CAPACITANCE (a, b, & c) 

TIME STEP SIZE (a&c) 

CATHODE DIAMETER (a, b, c & d) 
ANODE DIAMETER (a, b, c, & d) 
CATHODE HEIGHT (b, c, & d) 

ANODE HEIGHT (b, c, & d) 

GAS PRESSURE (a, b, c, & d) 


T 



REQUIREMENTS FOR A HIGH POWER SWITCH: 


FAST RISE TIME 

HIGH CURRENT CAPABILITY 

HIGH REPETITION RATE 

LOW WEAR OF ELECTRODE & INSULATORS 

(LONG LIFE) 


SWITCH COMPUTER 
SIMULATION MODEL 

(2-D, TRANSIENT 
MODEL) 



OUTPUT PARAMETERS: 

• PLASMA CURRENT (a, b, & d) 

© PLASMA DENSITY PROFILE (b. & d) 

• GAS DENSITY PROFILE (a & d) 

• PLASMA AXIAL VELOCITY PROFILE 
(a & c) 

• PLASMA RADIAL VELOCITY PROFILE 
(a & C) 

• PLASMA ION TEMPERATURE PROFILE 
(a, b & d) 

• PLASMA ELECTRON TEMPERATURE 
PROFILE (a, b & d) 

• MAGNETIC FIELD (a, b & c) 


PROTOTYPE SWITCH 
DESIGN 


*LETTERS IN PARENTHESIS DENOTE THE SWITCH 
REQUIREMENTS WHICH CAN BE DETERMINED BY DIRECT 
OR INDIRECT CONTRIBUTIONS OF PARAMETERS. 


ho 

ho 


Figure 5. Simulation Procedure for a High Power Inverse-Pinch Switch 





Figure 6. Switch Computer Simulation Model -Flow Chart. 



DEFINE: MODEL CONSTANTS 
AND INPUT PARAMETERS 



RESET: NEW INPUT 


PARAMETERS 


SET: TIME N=l, NMAX 
INITIAL CONDITIONS 


SET: DIRECTIONS 
I (RADIAL) = 1 
J (AXIAL) = 1 


GENERATE : TRI-DIAGONAL MATRIX 

ELEMENTS A , B , C , & D FOR 
s s s s 


• DENSITY EQ. 

• RADIAL VELOSITY EQ. 

• AXIAL VELOSITY EQ. 

•ION TEMPERATURE EQ. 

• ELECTRON TEMPERATURE EQ. 

• MAGNETIC FIELD EQ. 

• ELECTRIC CIRCUIT EQ. 


CALL: SUBROUTINES 

• VISCOUS PARAMETERS 

• CONDUCTION PARAMETERS 

• RESISTIVITY 

-• PLASMA DENSITY 

• COULOMB LOGARITHM 

• COLLISION RELAXATION TIME 

• BREMSSTRAHLUNG 

• MAGNETIC VISCOSITY . 


Y 


CALL: SUBROUTINE TRI-DIAGONAL 


SOLVE: MODEL EQUATIONS 


GENERATE: DATABASE 



" « N0 <r^ 


FOR PROTOTYPE DESIGN _*/ 

1 = 1 + 1 

I > IMAX 

> ^ 


j 

^ 



J = J + 1 


J > JMAX 


''OUTPUT PARAMETERS'^ no 
1EET DESIGN CRITERIA/” 


t YES 

SWITCH: DIRECTIONS 

I (AXIAL) - J 
J (RADIAL) = I 


N > NMAX 


23 







OBSERVATION WINDOW 



ANODE 

TERMINAL 


Figure 7. Cross-Section of a New Closing Inverse-Pinch 
Switch. 
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Figure 8. Cross-section of new Closing/Opening switch. The bottom 
portion of the C/0 switch consist of the inner and outer 
electrodes and HCP plasma-puff trigger mechanism for the 
closing phase. The top portion of the C/0 switch consist 
of the inner and outer electrodes and magnet for the pinch 
of plasma current. The opening switch uses m=0 instability 
to cut off the plasma current. 


LENGTH OF ANODE, 


TIME = .1 fx sec 

TIME = .8 ^4sec 

TIME = 3.1 Msec 



RADIUS 


Figure 9. Current sheet location along the radius at three- selected 
times. The number 1 (r-^ = 0.0365 m) corresponds to the l 
inner electrode (anode) and 13 (r Q = 0.0643 m) corresponds 
to the outer electrode (cathode) . The gas pressure was 
0.7 mTorr and the system inductance was 100 nH. The current 
sheet at 0.1 usee almost remains on the initial location. 


LENGTH OF ANODE, 


TIME = 
TIME = 


.1 /xsec 
.8 Msec 
3.1 Msec 



RADIUS 


Figure 10. Current sheet locations at 3 selected times while the 
gas pressure was maintained at 9 mTorr and the system 
inductanct at 100 nH. The numbers 1 and 13 are the 
inner and outer electrodes. The current sheet at 0.1 
y sec still stays on the initial postion. 
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SWITCH CURRENT, 
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AXIAL VELOCITY, M/SEC 


— RADIUS = 2 

RADIUS - 4~ : . 

— RADIUS = 6 

RADIUS = 8 

RADIUS = 10 



TIME, /XSEC 


i ■ • 

i • 

i Figure 14. Axial velocities of current sheet at 5 different radial 

| positions during the discharge period. The gas pressure 

| was 0.7 mTorr and the system inductance was. 100 nH. 
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TIME, MSEC 


Figure 16. Axial velocities at 5 selected radial positions during 
the discharge period, p = 0.7 mTorr and L = 700 nH 




TIME. yuSEC 


Figure 17. Radial velocities of current sheet at 5 different 
radial positions. The gas pressure was 0.7 mTorr 
and the system inductance was 100 nH. 









R 


Figure 20, Axial velocities along the i 
discharge times. The gas pi 
system inductance was 100 nf 
Radius axis are the inner ( 
(r G = 0.0643 m) electrodes 
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AXIAL VELOCITY, M/SEC 



Figure 23. Axial velocities at 3 selected pressures and at 

0.8 psec after the discharge began. L = 100 nH. 

sys 
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AXIAL VELOCITY, M/SEC 





RADIAL VELOCITY, M/SEC 



RADIUS 


Figure 26. Radial velocities at 3 selected times along the radial 
direction, p = 8 mTorr and L ■«* 100 oH. 





LEAKAGE CURRENT, 


DENR 

DENR 

DENR 


.7073999E-03 TORR 
•4337400E-02 TORR 
•7967400E-02 TORR 



Figure 27. Leakage current during the discharge period at 3 selected 
pressures, while L s „ s = .100 nH. It shows no effects on 
leakage current due to the variation of gas pressure. 


LE = .1 000000E-06 

LE = .4000000E-06 

LE = .7000000E-06 



Figure 28. Leakage currents at selected system inductances 
p = 0 . 7 mTorr . 
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DISCHARGE VOLTAGE, 


LE 

LE 

LE 


.1 OOOOOOE-O 6 
.4000000E-06 
.7000000E-06 



10.00 1 L 

0 .5 




J_J I 


1.0 1.5 


2.0 2.5 3.0 

TIME. /jSEC 


3.5 4.0 


Figure 29. Discharge voltage variations at 3 selected system 
inductances during the discharge period. 


DENR «= .7073999E-03 TORR 

DENR « .1 91 7400E-02 TORR 

— DENR = .31 27400E-02 TORR 
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DENR « .7967400E-02 TORR 

DENR = .9177400E-02 TORR 



TIME, yuSEC 


Figure 30 • Switch inductances at 8 selected pressures during the 
discharge period. L S y g = 100 nH. 
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APPENDIX A 


A new type of an inverse switch was designed and built during the period 
of the contract. The switch was installed onto a high power source in the 
laboratory of Space Technology Branch- Its trigger mechanism (HCP plasma 
puff) and high power load are currently being tested to determine the uniformity 
of trigger puff and voltage hold-off. 

The drawings of the switch parts are attached with an assembly drawing. 

As a result of testing, minor changes or corrections of the switch configura- 
tion were made for possible imporvement of switch performance. Therefore, 
assembly and parts drawings attached here might show some differences because 
the recent changes were not included in the drawings. 

Figure A-l shows the assembly drawings of a new switch. The numbers 
on the drawings designates the parts. The list of drawings follows: 

Fig. A Assembly Drawing 

Fig. A-l Observation View Window 

Fig. A-2 Anode Electrode Head-disk 

Fig. A-3 Switch Ghamber Housing (top) 

Fig. A-4 HCP Plasma Trigger Insulator 
Fig. A-5 Switch Ghamber Housing (bottom) 

Fig. A-6 Cathode Electrode 

Fig. A-7 Switch Support Substrate 

Fig. A-8 Switch Insulator Between Anode and Cathode 

Fig. A-9 Anode Electrode Column 

Fig. A-10 Anode Electrode Cable Connector 

Fig. A- 11 Electric Coaxial Cable Connector Frame 
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Fig. 

Fig. 


A- 12 Assembly of HCP Plasma Puff Trigger 
A-13 Connectors between View Window (A-l) and substrate 
(A-7). (not shown in Fig. A). 
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Figure A. Assembly drawing of a New Inverse-Pinch Switch. The numbers in the 
circles show part serial number which is the same as the figure 
number 

Scale: 1/1 . 
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Figure. A-5. Switch Chamber Housing: 
Material : Plexiglas 

Scale: 1/2 


Bottom part 



i 


Figure. A- 


o 

=P. 
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6. Cathode Electrode 
Material: Brass 

Scale: 1/2. 
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Figure. A-8. Switch Insulator 
Material : Macor 

Scale: 1/2 
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